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Abstract 



We recalculate the pion electromagnetic form factor based on the pertur- 
bative QCD formalism that includes the Sudakov resummation. We take into 
account the evolution of the pion wave function in b, which represents the 
transverse extent of the pion. An infrared enhancement is observed when 1/6 
approaches Aqcd- We propose to freeze the evolution of the wave function at 
some scale above Aqcd, such that perturbative contributions are stablized. 
Our predictions are consistent with experimental data, and insensitive to the 
variation of relevant parameters. 
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I. Introduction 



The applicability of perturbative QCD (PQCD) to exclusive processes 
has been a controversy for almost two decades |l|, |^. Although there is a 
general agreement that PQCD can successfully make predictions for exclusive 
reactions as momentum transfers go to infinity, it remains unclear whether 
experimentally accessible energy scales are large enough to justify these pre- 
dictions. A progress was made recently by Li and Sterman who proposed 
a modified PQCD formalism with Sudakov suppression taken into account. 
They found that nonperturbative contributions are suppressed significantly 
by the Sudakov effects, and the applicability of PQCD can be extended down 
to few GeV scales. However, their analysis of the pion electromagnetic form 
factor did not concern the match with experimental data, such that a con- 
vincing piece of justification of their approach is missing. In this paper we 
shall further include the evolution of the pion wave function into the modified 
PQCD formalism, and demonstrate that our predictions explain the data. 

The PQCD theory for exclusive processes was first proposed by Brodsky 
and Lepage [^f. They argued that an exclusive process, such as the pion 
form factor, can be factorized into two types of subprocesses: wave functions 
which carry the nonperturbative information of the initial and final state pi- 
ous, and a hard amplitude which describes the scattering of a valence quark 
of the pion off the energetic photon. The former is not calculable in per- 
turbation theory, and needs to be parametrized by a model, to be derived 
by nonperturbative methods, such as QCD sum rules and lattice gauge the- 
ory, or to be determined by experimental data. The latter, characterized 
by a large momentum flow, is calculable in perturbation theory. Combin- 
ing these two types of subprocesses, predictions for the pion form factor at 
experimentally accessible energies scales can be made. 

According to this naive picture, the factorization formula for the pion 
form factor graphically represented by Fig. 1(a), is written as 



= —2Pi ■ P2 is the momentum transfer from the photon. Pi (P2) being 
the mometum of the incoming (outgoing) pion. /i is a renormalization and 
factorization scale, below which QCD dynamics is regarded as being non- 
perturbative and absorbed into the pion wave function 0, and above which 




(1) 
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QCD dynamics is regarded as being perturbative and absorbed into the hard 
amphtude H. /x) gives the probabihty of a valence quark carrying a 
fractional momentum xP in the parton model at the scale /z. Usually, the ix 
dependence of is neglected, i.e., /i) = H is obtained by comput- 

ing quark-photon scattering diagrams. 

To make predictions for F^r, one substitutes the lowest-order expression 
of H derived from Fig. 1(b), 



H{xi,X2,Q,fJ.) 



XiX2Q'^ 



(2) 



and the asymptotic pion wave function 



kAS, 



X) 



3/. 



:X(1 — X) 



(3) 



into Eq. (P, where ^^(/i) is the running coupling constant, Cp = 4/3 the 
color factor, Nc = 3 the number of colors, and /^r = 0.133 GeV the pion decay 
constant. Because (j)^^ peaks at x = 1/2, the main contributions to Eq. (|^) 
come from the region with intermediate x. Higher-order corrections to H 
then produce logarithms of the type ln((5^//i^), which may be so large as 
to spoil the perturbative expansion. To eliminate the logarithms, a natural 
choice of /i is /i = Q. It is easy to find that Q^F^(Q^) ~ 0.12 GeV^ for ~ 4- 
10 GeV^ is only 1/3 of the data ~ 0.35 GeV^ [|, 0, |, ||. This contradiction 
implies that the pion form factor has not yet become completely asymptotic 
at experimentally accessible energies. 

Hence, one may resort to a preasymptotic model, the Chernyak and Zhit- 
nitsky (CZ) wave function derived from QCD sum rules [p!0|] . 



15/. 
V2N, 



x(l-x)(l-2x)^ 



(4) 



which possesses maxima at the end points x ^ and x —>■ 1. The CZ wave 
function enhances the predictions for Q'^F^^^Q'^) to 0.30 GeV^, and improves 
the match with the data. Howerer, such a success was criticized by Isgur 
and Llewellyn Smith |TT| and by Radyushkin []T2|: Since 0*-^^ empha- 



and by Radyushkin []T2|: Sin( 
sizes the contributions from small x, higher-order corrections to H in fact 
produce logarithms like \n{xiX2Q'^ / ^j^'^). Hence, the natural choice of is 
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/i = y/xiX2Q, which is consistent with jj. = Q employed above for intermedi- 
ate x, but results in a large coupling constant as{^/x^X2Q) > 1 at small x, 
indicating that the end-point regions are nonperturbative regions. A care- 
ful analysis revealed that the enhancement of the predictions for Q^F^^IQ'^) 
is due to the amplification of the nonperturbatve end-point contributions. 
As a consequence, the perturbative calculation loses its self-consistency as a 
weak-coupling expansion. 

The above discussions hint that the end-point regions should be treated 
in a different way. In general, a valence quark in the pion can carry a small 
amount of transverse momenta kx, which then flow through the hard scat- 
tering. Therefore, the hard gluon propagator is written as 



1 1/4 



XiX2Q'^ + XiX2Q'^ \ XiX2Q'^ 



(5) 



where the first term corresponds exactly to the asymptotic H in Eq. @, and 
the second is suppressed by powers of at intermediate x. If the end-point 
regions are not important, one may drop the second term, and integrate 
out the kT dependence in the wave functions, arriving at the leading-power 
factorization formula in Eq. (|l]). However, the end-point difficulities indicate 
that the higher-power effects are crucial, and should be kept at the outset as 
deriving the perturbative expression for the pion form factor. 

The introduction of kx leads to three observations immediately. First, 
the pion form factor becomes a two-scale {Q and kx) problem, and the re- 
summation technique is required to organize the large logarithms ln{Q^/k^) 
from radiative corrections to the wave functions. Second, the process must 
be analyzed in the Fourier transform space of k^, denoted by b which 
is regarded as the transverse separation between the valence quarks of the 
pion. Third, the quantity 1/6, now considered as one of the characteristic 
scales of the hard amplitude, should be substituted for the argument of as if 
1/6 > ^/xiX2Q. 

Li and Sterman found that the resummation of the large logarithms sup- 
presses the elastic scattering at large spatial separation. This property, called 



Sudakov suppression [|13| , p!4[ makes nonperturbative contributions from 
large b, no matter what x are, less important, and improves the applicability 
of PQCD at few GeV scales. However, their discussion was restricted to 
the PQCD apphcability, such that the results of the pion form factor do not 



4 



match the data. In the present work we shall aim to explain the experimen- 
tal data by including a new piece of information, the evolution of the wave 
function, into the modified PQCD formalism. Our predictions are found to 
be dominated by perturbative contributions for > 4 GeV^, and in good 
agreement with the data. 

A brief review of the standard asymptotic expression for the form factors 
[^, [l^, |16|, 0, |l9l may help to motivate our viewpoint. The factorization 



formula for the pion form factor including the evolution of the wave function 
is derived in Sect. II. Numerical results are given in Sect. III. Section IV is 
the conclusion. 



II. The Factorization Formula 

In this section we review factorization theorems including Sudakov ef- 
fects for the pion form factor, which is expressed as the convolution of a hard 
amplitude with wave functions. We investigate radiative corrections to the 
quark-photon scattering diagrams shown in Fig. 2, and explain how they are 
factorized into the convolution factors. There are two types of important 
corrections: collinear, when the loop momentum is parallel to the incoming 
or outgoing pion momentum, and soft, when the loop momentum is much 
smaller than the momentum transfer . Each type of the important cor- 
rections gives rise to large single logarithms. They may overlap to produce 
double logarithms in some cases. In axial gauge the two-particle reducible 
diagrams, like Figs. 2(a) and 2(b), give double logarithms from the overlap of 
collinear and soft divergences, while the two-particle irreducible corrections, 
like Figs. 2(c) and 2(d), give only single soft logarithms. 

It can be shown that soft divergences cancel between Figs. 2(a) and 2(b), 
as well as between 2(c) and 2(d), in the asymptotic region with 6 — 0. There- 
fore, reducible corrections are dominated by collinear divergences, and can 
be absorbed into the pion wave funtion P, which involves similar dynamics. 
Irreducible corrections, due to the cancellation of soft divergences, are then 
absorbed into the hard scattering amplitude H . Based on the above rea- 
soning, the factorization formula for the pion form factor is written in the h 
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space as 



FAQ' 



dxidx2 



X H{xi, X2, b, Q, fi)V{xi, b, Pi, /i) 



(6) 



where H is the Fourier transform of H. 

We choose the Breit frame such that Pi' = P2 = <5/"\/2 and all other 
components of P's vanish. After resumming the double logarithms, we obtain 



V{x, 6, P, /i) = exp 6, Q) - s(l - x, 6, Q)] V{x, 6, /i) 



The exponent b,Q), ^ = x and 1 — x, is expressed as [|14 



s{^,b,Q) 



eQ/^/2 dp 
1/b p 



In 



^V2p, 



A{as{p)) + B{as{p)) 



(7) 



(8) 



where the anomalous dimensions A to two loops and B to one loop are given 
by 



A 
B 



Cf— + 



67 TT^ 
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2a, 

3 TT 



In 
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27 



+ In — 



(9) 



7£; being the Euler constant. The two-loop running coupling constant, 

_ 4 _ 16/3i lnln(/iVA^) 
~ /3oln(/i7A2) ~ /33 ln2(/i7A2) ' 

with the coefficients 



/3o 



33 - 2n 



/ 



153 - 19nf 
6 



(10) 



(11) 



and the QCD scale A = Aqcd, will be substituted into Eq. (||). We require 
the relation of the involved scales ^Q/a/2 > 1/6 > A as indicated by the 
bounds of the variable p in Eq. QCD dynamics below 1/b is regarded as 
being nonperturbative, and absorbed into the initial condition V{x,b,^). 
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The functions V and H still contain single logarithms from ultraviolet 
divergences, which can be summed using the RG methods [14|. The large-6 
behavior of V is then summarized as 

6, P, /i) = exp ' ''' - 1- f - 



-s{x, b, Q) - s{l -x,b,Q) -2 f — 7g(a^(/i)) 

Jl/b u 



l/b fi 

xV{x,b,l/b), (12) 

with 7g = —as/n the quark anomalous dimension in axial gauge. To con- 
centrate on the evolution of V in l/b, we ignore the intrinsic b dependence 
denoted by the argument b ||2^, that is, we adopt V{x, b, l/b) = 0(a;, l/b). 

Since the pion form factor, as a physical observable, is /i independent, the 
RG analysis apphed to H gives 



H{xi,b,Q, jj,) = exp 



— 7q(a,(/i)) 



(13) 
(14) 



xH{xi,b,Q,t) , 

where t is the largest mass scale involved in the hard scattering, 

t = raax{^/xlX2Q, l/b) . 

The scale yJx^x^Q is associated with the longitudinal momentum of the hard 
gluon and l/b with the transverse momentum. At lowest order, H is given 
by 

l^nasCpXYQ^ 



[xiQ^ + kfy)(xiX2Q^ + (kir - \^2t)^) 
l^'KasCp 

x.x^Q' + (kiT - k2T)2 ' ^ ^ 

where we have neglected the transverse momentum in the numerator. In the 
second form, we have further neglected the transverse momentum associated 
with the virtual fermion lines in H, which are linear rather than quadratic 
in X. Because Eq. (|T^) depends on the combination of transverse momenta, 
H involves only a single b, which has been made explicit in Eq. (^. 

Combining all the above ingredients, we arrive at the factorization for- 
mula for the pion form factor, 

F^(Q2) = levrCj 



dxidx2 / bdb(j){xi,l/b)(p{x2,l/b)as{t)Ko{y/x^Qb) 
Jo 

(17) 



X exp[-S{xi,X2,b,Q)] . 
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where the complete Sudakov exponent is given by 



S{xi,X2,b,Q) 



2 f 
J2[six.,,b,Q) + sil-x,,b,Q)]+4: / 

i=i ^1 



/b Ji 



%{asm. (18) 



Kq is the modified Bessel function of order zero, which comes from the Fourier 
transform to h space of the gluon propagator. Note that we insert the two- 
loop as in Eq. ([T0|) into the above integral of 7g, while it is the one- loop as 
that was employed in The exponential decreases rapidly at large 6, 
leading to Sudakov suppression mentioned before. The essential advantage 
of Eq. (p!7D, compared to Eq. (|1|), is then the extra h dependence in the hard 
scattering. If h is small, radiative corrections with the argument of as set to 
t will be small, regardless of the values of x. Of course, when h is large and 
XiX2Q'^ is small, radiative corrections are still large. However, we shall show 
that e""^ suppresses the contributions from this region, and the perturbation 
theory becomes relatively self-consistent. 

In Eq. ([17| ) the evolution of with 1/6 is written as ||2l| 



1/6) 



v/2Ar 



■Ml 



X] 



(^nC^ (2x — 1 J 



Ln=0 



7n' 



(19) 



where the coefficients a„ and the exponents 7„ are 

5(2n + 3) ^ 



(n + 2)(n + 1) 
4 



7n 



33 - In 



f 



n+l 1 



k=2 



k {n + l){n + 2) 



with Jo = 2/15, Ji = 0, /2 = 8/35 and 



Gagenbauer polynomials are C( 
Then Eq. (|T^ reduces to 

3 



3/2 



X) 



or n > 3. 

(x) 



(20) 
(21) 

The relevant 



1 and C2^\x) = (15/2) (x2 - 1/5). 



0(x, 1/6) 



V2iV, 



/^x(l - x) 



1 + [5(1 - 2xf - 1] 



50 
81 



where we have changed the argument of as from 1/6 to 

p, = max(l/6, m) , 



(22) 



(23) 
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with m a mass scale greater than the initial value no = 0.5 GeV. This modifi- 
cation, implying that 1/6) is "frozen" at (f){x,m) as 1/6 < m, guarantees 
Eq. to evolve toward the correct direction. Without this freezing scale, 
1/6) is ill-defined for 1/6 — A, and large infrared contributions to the 
pion form factor will be observed, m will be treated as a parameter of order 
/io, and determined by the data fitting. Obviously, Eq. ( ^2]) approaches the 
asymptotic wave function (p^^ in Eq. (|]) as 1/6 ^ cxd, and the preasymptotic 
wave function 0*"^ in Eq. (|) as 1/6 — > /iQ. We then expect that the pion 
form factor derived from Eq. (p2D is located at between those from (p^^ and 



lCZ 



We are also interested in the full expression of H in Eq. (0), except 
for the approximate one in Eq. (JT^). Without neglecting the virtual quark 
transverse momemta, the parallel formula to Eq. ([T7|) involves a double-6 
integral: 

rl poo 

Ftt{Q^) = IQttCf / dxidx2 / bidbib2db2(j){xi, l/bi)(f){x2, 1/62) 
Jo Jo 

X as{t)Ko{y/xiX2Qbi) exp [-^(xi, X2, 61, 62, Q)] 



X [^(61 - 62)iro(^Q6i)Jo(v^Q62) + 

^(62 - b,)Koi^,Qb2)IoiV^iQbi)] , (24) 

with Iq the modified Bessel function of order zero, and the complete Sudakov 
exponent 

S = Vs(xi,6i,(5) + s(l -Xi,6i,Q) + 2 / — 7q(a^(/i)) , (25) 
and the hard scale 

t = max [^/ila^Q, 1/61, 1/62] . (26) 



A freezing scale is not necessary for Eq. (p^, since the extra degrees of free- 
dom 62 moderate the infrared enhancement from the low end of the evolution 
of the wave functions. This will be explicitly demonstrated in next Section. 



III. Numerical Results 
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In the numerical analysis the exponential e""^ will be treated as follows. 
In the small b region higher-order corrections should be absorbed into the 
hard amplitude, instead of into the wave functions, giving their evolution in 
b. Hence, we set any Sudakov factor e~'^(^'^''5) to unity for < 1/6. 

In this vein, we also set to unity, whenever it exceeds unity in the small 6 
region. As 6 increases, decreases, reaching zero at 6 = 1/A. Suppression 
in the large 6 region is weaker for smaller Q, where Sudakov effects are mild. 
The above typical behavior of has been shown in ||^. It is expected that 
the suppression improves the self-consistency of the perturbative calculation 
of the pion form factor. 

We evaluate Q'^F.^{Q'^) employing the pion wave function (p with evolution 
in Eq. (^. We find, at = g GeV^, Q^F^{Q^) = 0.33 GeV^ for the freezing 
scale m = /io and Q'^F.„{Q'^) = 0.28 GeV^ for m = 2/io, the latter being only 
15% smaller than the former. It implies that our formalism is insensitive to 
the choice of m. If a freezing scale was not introduced, the predictions are 
almost enlarged by a factor 6, which is attributed to the infrared enhancement 
from the region with 1/6 < /zq, i-e., to the misuse of Eq. (^2]). We adopt the 
natural choice m = hq in the analysis below. 

As to the self-consistency of the perturbative calculation, we propose to 
analyze Eq. (0), not by cutting off the x integrals near their endpoints 



but by testing the sensitivity of the integral to the large 6 region. We expect 
that, because Eq. (|lD is the true asymptotic behavior, as Q increases, the 6 
integral will become more and more dominated by the small 6 contributions, 
for which the effective coupling is small. To see how the contributions to 
Eq. (|T^ are distributed in 6 space under Sudakov suppression, the integration 
is done with a variable cut off in 6, be- Typical numerical results are displayed 
in Fig. 3 for the use of in Eq. (|22|) with A = 0.15 GeV. The curves, 
showing the dependence of on 6c, rise from zero at be = 0, and 

reach their full height at 6^ = 1/A, beyond which we consider any remaining 
contributions as being truly nonperturbative. As anticipated, we observe a 
faster rise as Q increases. 

To be quantitive, we may consider the cutoff up to which half of the 
whole contribution has been accumulated. For = 4 and 9 GeV^, 50% of 
Q'^F^{Q^) comes from the regions with 6 < 4.0 GeV"^ [a^(l/6)/7r < 0.44] 
and with 6 < 3.2 GeV~^ [as^l/b)/^ < 0.37], respectively. A liberal standard 
to judge the relevance of the perturbative method would be that 50% of the 
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results come from the region where as(l/6)/'7r is no larger than, say, 0.5. 
From this point of view, PQCD begins to be self-consistent at ~ 4 GeV^. 
This conclusion is consistent with that drawn in where the CZ wave 
function 0*^^ and the asymptotic wave function cf)"^^ were employed. 

Results of Q^F^iQ^) for the use of 0^^, and in Eq. (H), as well 
as the experimental data ||^, 0, H, are shown in Fig. 4. As expected, the 
curve corresponding to Eq. (^21) is located at between the other two. It is 
known that the pion form factor behaves like asymptotically. The data 
points of Q^Ft.^Q'^), arising from low and reaching a pleatau at ~ 2 
GeV^, reflect this asymptotic behavior. Hence, the point at ~ 10 GeV^, 
jumping to a value about twice of the pleatau, seems not plausible, and is 
not considered here. In this sense, our predictions are in a good agreement 
with the data for > 4 GeV^. 

For < 4 GeV^, all the curves rise rapidly and deviate from the data, 
indicating the dominance of the nonperturbative contrubutions. Note that 
the predictions from (p^^ and 0*"^ are larger than those obtained in This 
is due to the insertion of the two- loop as into the integral of in Eq. (p!8|). 
To examine the stability of our predictions to the variation of the scales A 
and /io, we compute Q'^F^.^Q'^) for A = 0.1 and 0.2 GeV {fio = 0.5 GeV), and 
for /io = 0.4 and 0.6 GeV (A = 0.15 GeV) at = g GeV^. Their values 
are listed in Table I. Apparently, the difference from those corresponding to 
A = 0.15 GeV and /iq = 0.5 GeV is small. 

At last, we evaluate Q'^Ft^{Q^) using the double-6 factorization formula 
in Eq. (2i), and the results are presented in Fig. 5. Those derived from the 
single-6 formula in Eq. (|l^) and the data are also displayed for comparision. 
We find that the two curves are close to each other especially in the region 
with large ~ 5-10 GeV^. The consistency of the predictions from Eq. (p^) 
with the data is nontrivial, since it does not involve an adjustable freezing 
scale m. The above analyses indicate that our scheme of truncating the in- 
frared enhancement from the evolution of the wave function indeed makes 
sense. In the future study of the proton form factor, we shall take the same 
approach in order to simply the analysis, that is, neglect the transverse mo- 
menta carried by the virtual quarks, and introduce a freezing scale for the 
proton wave function. 



IV. Conclusion 
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In this paper wc have performed a very thorough investigation of the 
pion form factor based on the modified PQCD formahsm. Both the two ap- 
proaches with the single-6 formula plus the introduction of a freezing scale 
and with the double-6 formula lead to results that match the experimental 
data. We have also confirmed that our predictions are dominated by per- 
turbativc contributions for > 4 GeV^, and insensitive to the variation of 
the parameters A and Hq. All these conclusions imply that our formalism 
is rehable for the study of exclusive QCD reactions at intermediate energy 
scales. We shall extend it to more complicated processes, such as the proton 
form factor and Compton scattering, which will be published elsewhere. 
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Table I. The values of Q^F^{Q^) for different A and /xq at = 9 



A(GeV) 


g^F,(Q^)(GeV^) 


/^o(GeV) 


g^F,(Q^)(GeV^) 


0.10 


0.2974 


0.4 


0.3114 


0.15 


0.3309 


0.5 


0.3309 


0.20 


0.3659 


0.6 


0.3333 
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Figure Captions 

Fig. 1. (a) Factorization of the pion form factor. Tlie symbol x represents 
the photon vertex, (b) Basic scattering diagrams. 

Fig. 2. Radiative corrections to the basic scattering diagrams. 

Fig. 3. Dependence of Q^F^{Q^) on be for the use of in Eq. (H) at = 4 
GeV^ (dashed hne) and at = 9 GeV^ (sohd hne). 

Fig. 4. Dependence of on Q"^ for the use of (p^^ (dotted line), of 

(j)'^^ (dashed line), and of (j) in Eq. (^) (solid line). 

Fig. 5. Dependence of Q'^F^^^Q^) on for the use of (f) in Eq. ( ^2]) derived 
from Eq. (|T^ (dashed line) and from Eq. (^) (solid line). 
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